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We introduce a new parameterization of four-fermion matrix elements which does not involve quark masses and 
thus allows a reduction of systematic uncertainties in physical amplitudes. As a result the apparent quadratic 
dependence of e'/e on m s (/i) is removed. To simplify the matching between lattice and continuum renormal- 
ization schemes, we express our results in terms of Renormalization Group Invariant B-parameters which are 
renormalization-scheme and scale independent. As an application of our proposal, matrix elements of AI — 3/2 
and SUSY AF — 2 (F = S, C, B) four-fermion operators have been computed. 



1. Introduction 

Since the original proposals of using lattice 
QCD to study hadronic weak decays sub- 
stantial theoretical and numerical progress has 
been made: the main theoretical aspects of the 
renormalization of composite four-fermion oper- 
ators are well understood ^ |(|; the calculation 
of K°-K° mixing, relevant to the prediction of 
the CP- violation parameter e, has reached a level 
of accuracy which is unpaired by any other ap- 
proach [Q; increasing precision has been gained 
in the determination of the electro- weak penguin 
amplitudes necessary to the prediction of the CP- 
violation parameter e'/e p|-|lO|]; finally matrix ele- 
ments of AS = 2 operators which are relevant to 
study FCNC effects in SUSY models have been 
computed |],|l(|. Methods and symbols used in 
this talk and all the results we report are fully 
described in P,|To|]. 

2. Matrix elements without quark masses 

The analysis of K° — K° mixing with the most 
general AS = 2 effective Hamiltonian requires the 
knowledge of the matrix elements (K°\Oi\K°) of 
the parity conserving parts of the following oper- 
ators 

Ox = .5^(1-75)^^(1-75)^, 
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2 = ^(1-75)^^(1-75)^, 

03 = ^(1-75)^^(1-75)^, (1) 

04 = s a (l-j 5 )d a ^(1 + 75)^, 

5 = ^(1-75)^(1+75)^. 

On the lattice, matrix elements of weak four- 
fermion operators are computed from first prin- 
ciples. But, following the common lore, they 
are usually given in terms of the so-called B- 
parameters which measure the deviation of their 
values from those obtained in the Vacuum Sat- 
uration Approximation (VSA). For operators in 
([!]), the B- parameters are usually defined as 



(2) 



Ci 
3 



M 2 K f K 



Bi(fi). 



(K \O^)\K°) , 

where Ci = -5, 1, 6, 2 for (i = 2, ... , 5). In 
(g), (i\~°|Oi|-Kr ) is parameterized in terms of 
well-known experimental quantities and Bi(fi) 
(B K (p) = Bi(p)). On the contrary, (K°\Oi\K°) 
[i = 2, . . . , 5) depend quadratically on the quark 
masses in (g), while they are expected to remain 
finite in the chiral limit and depend only linearly 
on the quark masses. Contrary to fx, Mk, etc., 
quark masses can not be directly measured by ex- 
periments and the present accuracy in their deter- 
mination is still rather poor. Therefore, whereas 
for 0\ we introduce Bk as an alias of the matrix 
element, by using (||) we replace each of the SUSY 
matrix elements with 2 unknown quantities, i.e. 
the £?-parameter and m s +nid- To overcome these 
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Table 1 

Matrix elements in GeV A at the renormalization 
scale /i = 2 GeV in the RI scheme. In the first 
two columns the results obtained with the new pa- 
rameterization are given, while in the last two 
columns we show the results obtained in ref. ||/ 
with the "conventional" parameterization on the 
same set of data. 





New 


Old 




(3 = 6.0 
this work 


[3 = 6.2 
this work 


(3 = 6.0 
§ 


{3 = 6.2 
§ 


(Oi) 
Bi 


0.012(2) 
0.70(15) 


0.011(3) 
0.68(21) 


0.012(2) 
0.70(15) 


0.011(3) 
0.68(21) 


(O2) 
B 2 


-0.079(10) 
0.72(9) 


-0.074(8) 
0.67(7) 


-0.073(15) 
0.66(3) 


-0.073(15) 
0.66(4) 


(0 3 ) 
B 3 


0.027(2) 
1.21(10) 


0.021(3) 
0.95(15) 


0.025(5) 
1.12(7) 


0.022(5) 
0.98(12) 


(0 4 ) 
Bi 


0.151(7) 
1.15(5) 


0.133(12) 
1.00(9) 


0.139(28) 
1.05(3) 


0.133(28) 
1.01(6) 


(0 5 ) 
B 5 


0.039(3) 
0.88(6) 


0.029(5) 
0.66(11) 


0.035(7) 
0.79(6) 


0.029(7) 
0.67(10) 


<o 7 3/2 } 

B 3 7 /2 


0.019(2) 
0.65(5) 


0.011(3) 
0.38(11) 


0.020(5) 
0.68(7) 


0.014(5) 
0.46(13) 


(Of) 


0.082(4) 
0.92(5) 


0.068(8) 
0.77(9) 


0.092(19) 
1.04(4) 


0.087(19) 
0.98(8) 



problems, we propose the following new parame- 
terization of AS = 2 operators 



Mlf K B^X 



(3) 



(K \Oi^)\K ) = Ml.flBiin). 



The i?i(/i) parameters arc still dimcnsionless 
quantities and can be computed on the lattice 
by studying appropriate ratios of three- and two 
point functions pOj. By simply using them, we 



Table 2 

RGI Matrix elements 
Eq. (pi ) with a" /_4 . 



GeV 4 computed as 





(3 = 6.0 


(3 = 6.2 


{0 uoi } 


0.017(3) 


0.016(4) 


<0^> 


-0.051(7) 


-0.048(6) 




0.005(7) 


-0.004(7) 




0.072(3) 


0.063(6) 


(O^i) 


0.043(3) 


0.032(5) 



have eliminated any fictitious reference to the 
quark masses, hence reducing the systematic er- 
rors on the corresponding physical amplitudes. 
An alternative parameterization which has not 
been used in our numerical analysis, but may be 
very useful in the future, can be found in ]io| . 
The VSA and B-parameters are also used for ma- 
trix elements of operators which enter the AS = 1 
effective Hamiltonian. Notice that this "conven- 
tional" parameterization is the only responsible 
for the apparent quadratic dependence of e'/e on 
the quark masses. This introduces a redundant 
source of systematic error which can be avoided 
by parameterizing the matrix elements in terms of 
measured experimental quantities and therefore 
a better determination of the strange quark mass 
m s (/x) will not improve our theoretical knowledge 
of e'/e. In this work we have computed the ma- 
trix elements (it\0^ 2 \K) of the four fermion op- 
erators Z J 2 (i — 7, 8, 9) which contribute to the 
AI = 3/2 sector of e'/e. In fact in the chiral limit 

(7T7r|0^ 2 \K) can be obtained, using soft pion the- 



orems, from {n + \0'. 
masses, m s — md 
we find 

lim>+|0f 2 |if+) 



.3/2 



|if + ). For degenerate quark 
to, and in the chiral limit, 



lim (ir^ 

m — >0 

lim (ir^ 

m — >0 



-Mlfl lim BsOi) 
-Mlfl lim B a (jm) 



\0T\K + ) = 
\0 3 g /2 \K+) = \Mlfl lim BrOu) • 



In the limit m s = m,d complicated subtractions 
of lower dimensional operators are avoided for 
AI = 3/2 operators. This is not the case for 
AI = 1/2 operators which enter the determi- 
nation of e'/e : in this case the mixing with 
lower dimensional operators makes the computa- 
tion much more involved. A reliable lattice esti- 
mate of these matrix elements is still missing but 
encouraging preliminary results with domain- wall 
fermions have been presented in ref. 11 1. 



3. Renormalization Group Invariant Oper- 
ators 

Physical amplitudes can be written as 
(F\H eff \I) = (F\d(n)\I)-C(f,), (4) 
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where U(fx) = (Oi(/x), . . . , Ojsr(fi)) is the opera- 
tor basis (for example the basis defined in ([j]) for 
the AS = 2 effective Hamiltonian) and C(/i) the 
corresponding Wilson coefficients represented as 
a column vector. C(fi) is expressed in terms of 
its counter-part, computed at a large scale M, 
through the renormalization-group evolution ma- 
trix W\pL, M] 



where 



C{jj) = W[n,M]C(M), 



(5) 



where the initial conditions C(M), are obtained 
by perturbative matching of the full theory to the 
effective one at the scale M where all the heavy 
particles have been removed. W[fi, M) can be 
written as (see for example Jf2[) 



W[n,M] =M[/x]!7[/x,M]M- 1 [M] 



(6) 



where U = {a s {M) / a s {p))^o / 2 M is the 
leading-order evolution matrix and is a 

NLO matrix defined in |f2| which can be obtained 
by solving the Renormalization Group Equations 
(RGE) at the next-to-leading order. The Wilson 
coefficients C(n) and the renormalized operators 
O(fi) are usually defined in a given scheme, at a 
fixed renormalization scale /i, and they depend 
on the renormalization scheme and scale in such 
a way that only H e fj is scheme and scale inde- 
pendent. This is a source of confusion in the liter- 
ature, especially when (perturbative) coefficients 
and (non-perturbative) matrix elements are com- 
puted using different techniques, regularization, 
schemes and renormalization scales. To simplify 
the matching procedure, we propose a Renormal- 
ization Group Invariant (RG1) definition of Wil- 
son coefficients and composite operators which 
generalizes what is usually done for Bk and for 
quark masses IHI. We define 



(0)T ,„„ 

~7o /2ft> 



w [lA = M\fi] Mm)P° 
and using Eqs. @ and (0) we obtain 



(7) 



(8) 



The effective Hamiltonian 



K 



= o 



RGI 



■ c 



RGI 



can be written as 



(9) 



C 



RGI 



= w[M)C{M) , HG1 = 0(^)-w->] .(10) 



C RGI and O rgi are scheme and scale indepen- 
dent at the order we are working. Therefore the 
effective Hamiltonian is splitted in terms which 
are individually scheme and scale independent. 
This procedure is generalizable to any effective 
weak Hamiltonian. The B-parameters defined in 
eqs. (||) satisfy the same RGE as the correspond- 
ing operators and the RGI B-parameters can be 
defined as 



£f G/ = E^MM) 



-l 

ji 



(11) 
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Figure 1. Comparison of {K°\Oi\K°) in GeV 4 com- 
puted with the new (square) and the old (circle) pa- 
rameterization as a function of the strange quark 
mass used in to obtain the full matrix element 
from the B parameters. 



4. Numerical results 

All details concerning the extraction of matrix 
elements from correlation functions and the com- 
putation of the non-perturbative renormalization 
constants of lattice operators can be found in 
P,p|,|i"cfl . In this talk we report the results ob- 
tained in The simulations have been per- 
formed at p = 6.0 (460 configurations) and 6.2 
(200 configurations) with the tree-level Clover ac- 
tion, for several values of the quark masses and 
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for different meson momenta. The physical vol- 
ume is approximatively the same on the two lat- 
tices. Statistical errors have been estimated with 
the jacknife method. The main results we have 
obtained for AS — 2 and Al = 3/2 matrix el- 
ements and their comparison with the results in 
H are reported in Tables [l] and |[ In Figure [j] we 
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Figure 2. Values of (K°\Oi\K°) in GeV 4 computed 
at fi = 2 GeV in the RI scheme. Oi corresponds to 
the Standard Model AS = 2 operator. 

show the strong dependence of (K a \04\K°) on 
the strange quark mass when the "conventional" 
parameterization (||) is used, to be compared with 
the results obtained with the new parameteriza- 
tion. It is also evident that with the same set 
of data the new parameterization allows to deter- 
mine the matrix elements with smaller systematic 
uncertainties. Although we have data at two dif- 
ferent values of the lattice spacing, the statisti- 
cal errors, and the uncertainties in the extraction 
of the matrix elements, are too large to enable 
any extrapolation to the continuum limit a — > 
: within the precision of our results we cannot 
study the dependence of i?-parameters on a. For 
this reason, we estimate our best values of the 
£?-parameters by averaging the results obtained 
at the two values of /3 |lo]]. Since the results at 
j3 = 6.0 have smaller statistical errors but suf- 
fer from larger discretization effects, we do not 
weight the averages with the quoted statistical 
errors but take simply the sum of the two val- 
ues divided by two. As far as the errors are con- 



Table 3 

Matrix elements in GeV 1 at /i = 2 GeV in the RI 
scheme and their RGI values with a s f ~ . 



(Oi) 


RI 


RGI 


(Oi) 


0.012(3) 


0.017(4) 


(0 2 ) 


-0.077(10) 


0.050(7) 


(0s) 


0.024(3) 


0.001(7) 


(0 4 ) 


0.142(12) 


0.068(6) 


(0 5 ) 


0.034(5) 


0.038(5) 



cerned, we take the largest of the two statistical 
errors. Our best results are reported in Table || 
and are shown in fig. ^|. It is interesting to note, 
as expected from chiral perturbation theory, that 
matrix elements of AS = 2 SUSY operators are 
enhanced respect to the SM one ((Oi)) by a fac- 
tor 2 — 12 at /i = 2 GeV. Therefore, low energy 
QCD effects can enhance contributions beyond 
the Standard Model to ck [[l5||l6| l which, com- 
pared with the other SM predictions, becomes 
a promising observable to detect signals of new 
physics at low energy. The results for the analo- 
gous AC = 2 and AB = 2 matrix elements are 
reported in fl7| . 
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